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Abstract 

The topological vacuum structure of the two-dimensional CP n ~ l model 
for n = 3, 5, 7 is studied on the lattice. In particular we investigate the 
small-volume limit on the torus as well as on the sphere and compare with 
continuum results. For n > 5 , where lattice artifacts should be sup- 
pressed, the topological susceptibility shows unexpectedly strong deviations 
from asymptotic scaling. On the other hand there is an indication for a 
convergence to values obtained analytically within the limit n — > 00 . 
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1. Introduction 



The 2D CP n ^ 1 model [|l[ is a convenient play-model for testing non- 
perturbative methods aimed for application in asymptotically free field the- 
ories like QCD. Because of the existence of classical (multi-) instanton solu- 
tions it allows a semi-classical treatment @, • It can be considered within 
approximations based on the 1/n expansion 0, || as well as within the 
lattice approach [6 - 14]. 

One of the main quantities calculated within these different approaches 
is the so-called topological susceptibility \t , i- e. the zero momentum 
correlator of the topological charge density For quenched QCD the analogous 
quantity has an important physical meaning, since it is directly related to 
the mass of the rj meson (solution of the Ua(1) problem [fiTf ). 

However, in computing xt one meets serious difficulties both in the con- 
tinuum and on the lattice. In the continuum the semi-classical approach can 
be developed in a closed form only for the sectors of exact multi- (anti-) in- 
stanton amplitudes. The contribution of mixed instanton and antiinstanton 
configurations to the corresponding functional integral has to rely on model 
assumptions, like the dilute instanton gas or liquid models in QCD. Only in 
the small- volume limit, where the one-instanton contribution dominates, reli- 
able results [0 were obtained consistent with those of an 1/n expansion ||. 
Technically this achievement was possible by formulating the CP n ~ x model 
on a sphere S 2 . (The same holds for the pure SU(N C ) Yang-Mills theory 
treated on S A [jT7|| ). On the other hand, on the lattice there is the danger 
to encounter lattice artifacts, which were shown to spoil the continuum limit 
of the topological susceptibility at least for n < 3 J6|, [7|, |j. The analogous 
problem occurs in pure 577(2) lattice gauge theory when the geometric 
method to define the topological charge |T3 is directly applied to the Monte 



Carlo generated quantum fluctuations of the gauge field [19 



All these difficulties so far have prevented a real consistency check between 
the continuum approaches and lattice computations. Since new techniques 



for updating the complex vector fields are available [22[ we can study the 
lattice CP n_1 model for higher n , where lattice artifacts should be sup- 
pressed. Furthermore, we treat the model both on a torus and on a sphere. 
The latter case gives us the opportunity to investigate the topological struc- 
ture in the small-volume limit with the same boundary conditions as it was 
done in the continuum case flip]. 



In our previous papers [10 we concentrated on the CP 3 model. We 
used a lattice formulation of the model, where the Abelian gauge field in- 
volved was considered as a dynamical one. We found a situation which bears 
a great resemblance to that in 577(2) lattice gauge theory: With the geo- 
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metric method one finds an almost perfect asymptotic scaling behaviour for 
the topological susceptibility even on very small lattices, although there are 
dislocations expected to spoil the continuum limit. Thus, one is tempted 
to ask whether the scaling behaviour is accidental and could be disturbed if 
dislocations and other short range fluctuations are suppressed. 

Here we come back to another lattice formulation first considered by 
Petcher and Liischer |§ and apply to it the cluster updating algorithm 



12|j . The cases n = 3,5,7 will be treated. We shall investigate the 
topological properties and estimate the correlation lengths £ both for the 
two-dimensional torus T 2 and for the sphere S 2 . We are going to discuss 
the behaviour of the topological susceptibility down to volumes V ~ £ 2 
and to compare our lattice results with those relying on the 1/n expansion. 

Section 2 will collect all continuum results needed for our discussion. The 
lattice formulation on T 2 and S 2 , respectively, will be presented in Chapter 
3. Section 4 is devoted to the Monte Carlo results. The conclusions will be 
drawn in Section 5. 



2. The CP" 1 Model 

The CP™' 1 model on an arbitrarily curved two-dimensional manifold is 
defined by the following action [|TJ 

S= 2fJ dx ^^ D » Za DvZa ' (1) 

where g^ u denotes the metric tensor. The complex vector field z a (x), a = 
1, ...,n is subjected to the normalization condition 

z a z a = l . (2) 

stands for the total covariant derivative, i.e. it consists of the derivative 
in a curved space and in our case a U{\) gauge field A^, 

Dp — + iAp . (3) 

acts on an arbitrary vector in the curved space according to 

V„B V = d^B u - Y%B a , (4) 

where denote the Christoffel symbols. 

Expressing the gauge field A^ completely in terms of the z-field, 

A fl = -~(z a d tl z a -z a d»z a ), (5) 
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enables us to write 

s = Yj I d2x ^ 9 ' w ( d ^ + ) • ( 6 ) 

It is convenient to introduce projection matrices as follows || 

P a p{x) = z a (x)zp(x) (7) 

satisfying 

P = p\ p 2 = p, TrP = 1 Vz . (8) 
Expressing (P in terms of (0), we get 

5 = J fxy/ggT Tr [d,Pd u P] . (9) 

For our model the topological charge Q can be defined as 

Q = f d 2 xp(x) , (10) 



p{x) = — s^v V M A, 

Z7T 



where £12 = —£21 = 1 • This charge can be understood as the quantized 
"magnetic flux" through the surface. Expressing this "flux" in terms of the 
^-fields we get 

Q = -77- / d 2 X Sav df,z d v z . (11) 

Z7T J 

The following relation connects the action (P and the topological charge ([H]) 

77" 77 

S>j\Q\- (12) 

In case that (12) becomes an equality, we are led to the self-duality relation 

D"z = T^p£ua D a z . (13) 

For the 2-dimensional flat infinite space the (multi-) instanton solutions of 
(|13|) look as follows |6| 



z Pa[X) 

\p(x)\ ' 
Q 

Pa = C a Y[(xi - ix% - a 3 a ), C n = l. (14) 
3=1 



All the o? a are different from each other, and x\,Xi denote the coordinates 
in flat space. Note that the topological charge is determined by the zeros 
of the polynomials p a . In the case of our manifold being the sphere S 2 the 
coordinates X\,X2 have to be understood as the stereographic coordinates. 

The model is quantized using the functional integration method. Quan- 
tum expectation values of operators 0(z,z) have to be calculated as follows. 



(0(z,z)) 
DzDz 



Z~ l J DzDzO(z,z) exp(-S) 



2m 



5 ( Zp{x)zp{x) - 1 

vB=l 



(15) 
(16) 



with Z = J DzDz exp(— S) . The observables we are really interested in 
throughout this paper are the topological susceptibility \t and the corre- 
lation length £. The latter can be read off from the long-distance behaviour 
(cf. the next paragraph) of the correlation function 

2\ 



C(x,y) 



J2z a (x)z a (y) 



J2(z a (x)zp(y)) (z p (y)z a (x)} 

a,0 



1 

n 



The topological susceptibility is defined as follows. 

(Q 2 ) 



Xt= d x(p(x)p(0)) 



V 



(17) 



where V is the "volume" (i. e. the surface area) of the manifold under 
consideration (either the sphere S* 2 or the torus T 2 ). xt can be rewritten as 



Xt 



V 



E Q 2p Q 



(19) 



with Pq the probability for a field configuration to have the topological charge 
Q. On the sphere S 2 this formula can actually be exploited within the limit 
of a small volume where it gives 

2 



Xt 



V 



Pi 



(20) 



i.e., only the one-(anti-) instanton contribution turns out to be dominant. 
The result is 



Xt(V)A- M 2 s 



2d(n) 
VA 2 ms 



An 



n ln VA-ms 

71 4:71 



(21) 
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The coefficients C\(n) are numerically known, in particular 



Ci(3) 
C x (5) 
Ci(7) 



27.264 
3.963 
0.343 . 



(22) 



On the other hand, the model is 1/n-expandable [§J. The result within the 
leading order and for V = oo is 



Xt(oo) 



47m£ : 



+ — 



n- 



(23) 



where £ denotes the correlation length , which tends to a finite large n 
limit 



A/5 



4A 



q-tFQ) = Q.188 . 



(24) 



On the contrary, if the model is treated within the 1/n expansion on a sphere 
S 2 having a small finite surface area V , then Xt(V) behaves in accordance 
with the semi-classical result (2l~D ||. The limits n — > oo and V —>■ oo do 
not commute. 



It is worthwhile for our further considerations to replace Ams in (j2lD by 
^ obtained at V = oo 



2Ci(ti) 



V_ (01 

Ait 



2\ 2 



n , V /0.18 
— In — 

7T 47T V ^ 



(25) 



3. The Lattice Formulation 

First, let us describe the lattice constructions for the torus and the sphere. 
We decided to use simplicial lattices as the most convenient way for discretiz- 
ing curved manifolds. In the torus case it will allow us directly to compare 
with and, as we shall explain below, with 1/n expansion results. 

Torus T 2 : The plaquettes of the lattice are equilateral triangles with edge 
length a. A lattice point x is then described by 

x = xitii + x 2 n,2 
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where 

n x = a (1,0), (26) 

- - -(if) (27) 

The topology of the torus is established by imposing periodic boundary con- 
ditions on the functions defined on the lattice. 

/ {S+N 1 n 1 + N 2 (n 2 - = / (*) (28) 

where iVi and N 2 /2 are integers. The surface area V of the torus is 

V = -^-N^a 2 . 



Sphere S 2 : The sphere S 2 is discretized by a simplicial lattice in the 
following way. We start from a regular tetrahedron with corners placed at 
the surface of the sphere. Then we make several steps of refinement of this 
initial lattice. During each step we project the center of each given link onto 
the surface to introduce further sites. All sites, the old and new ones, are 
connected with their nearest neighbours by straight lines. After each step 
every simplex is replaced by four new simplices. After iV steps we have 
P = 2(4 N + 1) sites and L = 3(P — 2) links. Figs. 1, 2 show subsequent 
steps of refinement and the obtained link length distribution, respectively. 
In order to define a metric on these lattices, we consider the triangles to be 
flat. Identifying one of the corners i with the origin of the local coordinate 
system, we have for the triangle (i,j,k) |23j 



k 




the following metric tensor 



72/2 





where A is the area of the triangle and U are the lengths of the links 
corresponding to the unit vectors nj. 
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Lattice action: The most convenient way to discretize the lattice action 
uses its representation in terms of the projection matrices @. Thus we have 



77 

4 / A 

where the partial derivative has to be taken as follows: 

P(x + l^) - P(x) 



(30) 



I, 



M = l,2 



(31) 



can be rewritten as a sum over links (ij) between nearest neighbour 
sites i,j 

— n 

(32) 



with 

and the weight factors 



Sis^q (l - \z(i)z(j)\ 2 ) 



^kj^ki 



+ 



2A(i,j, k) 2A(i,j,k') 



(33) 



(34) 



The corresponding geometry is shown in the figure. 

J 




The action ( |3"2] ) is obviously U(l) gauge invariant. Note that in the case of 
the torus T 2 the weights are constant: 



T%3 ~ V3' 



(35) 



i.e., we obtain the action invented in 



Topological charge: For a given field configuration on the lattice, the 
topological charge is computed by using its definition according to 0. 

Q = X>(A) g(A)e(-ii), (36) 
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where 

/AXX Tr (P(k)P(j)P(i)) 

The sum runs over all (positively oriented) triangles. Equation ([37]) rewritten 
in terms of the z-field means 

/.„ / a \ \ z(i)z(J) z(J)z(k) z(k)z(i) 

exp(!2 ^ (A)) = mm\ mm\ mm\ m 

= racp(i(0(i,j)+0(j.*)+0(M))) (39) 

and finally 

27rg( A) = </>(i,j) + <j>(j, k) + <f>(k, i) + 2vrm( A) e (-vr, vr] , (40) 
with m(A) a properly chosen integer. 

Monte Carlo simulation: The expectation values ([15]) are computed nu- 
merically using the Monte Carlo method. We generate equilibrium field con- 
figurations by applying the one-cluster algorithm invented by U. Wolff 



This algorithm can be easily implemented for arbitrary n and works reason- 
ably fast, in spite of the fact that it does not remove critical slowing down 
T2|| . We choose an arbitrary lattice site Xq and collect adjacent lattice 



points into the cluster with the bond probability 

p = (1 - exp(-/5A^)) 6(AS tJ ) , (41) 
ASij = S tj (Rz(i), z(j) ) - S id ( z(i), z(j) ) , (42) 



with Sij the link action fl3~3|) and 

^ - { i w 

Having found all sites belonging to the cluster, the field vectors sitting at 
these sites are flipped. The flip operation R acting on the fields z is 
defined as follows: 

Rz = z - 2(rz)r, (44) 
where r is a random vector in the z-space, i.e., 



a=l 

Note that 



Yl r ^a = 1. (45) 
=1 

R 2 z = z , (46) 
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(Rz(i))(Rz(j)) = z(i)z(j). (47) 
Choosing the bond probability and the flip operation as above endows the 



algorithm with detained balance |22|, p35| . Fig. 3 shows the cluster point 
distribution for the CP 6 Model on the torus T 2 with (3 = 4.6 on a 100 x 100 
lattice. The average number of cluster points is 6516. 

In order to test the reliability of our program we compared the expecta- 
tion value of the link action (|33|) with analytical results. On the torus, the 
following weak coupling formula holds H] 



P » 1 {Sij) = - 



n 



3 V 2 P 
and for (3 = the strong coupling ones 

2 n- 1 



1 1 (// 



l)n 1 



(Sij) 



T 



T-- 



y/3 n 
- 1 



16^ I3 2 
on the torus T 2 , 



+ 



n 



L a 



on the sphere S 2 , where 
with L the number of links. 



(48) 

(49) 
(50) 



Fig. 4 shows the analytical curve fl48f) compared with our numerical values 
for Sij versus (3 for the CP 2 and CP A Model on the torus T 2 . For (3 = 0, we 
found 





(5«> 




Analytic 


Numeric 




((H) resp. ©) 




torus T 2 






CP^ Model 


.7698 


.7701 


CP 4 Model 


.9238 


.9231 


sphere S* 2 






r - = 1.5416 (N = 3) 






CP 2 Model 


1.0277 


1.0266 


CP 4 Model 


1.2333 


1.2318 



The purely statistical errors are omitted because they were found of order 
10~ 5 . 

In order to compute the correlation function, we used the improved esti- 
mator method []23|, |T2|, i.e., 

C{x) = {n + l)( ri e{x,c)e(0,c)Re[{z{x)z{0) 



(z(x)r)(rz(0)))(z(x)r)(rz(0)) 
10 



(51) 



where r denotes the reflexion vector ([0]) and 



„/ \ J 1 if x belongs to the cluster , . 

| otherwise. 

P and |c| are the number of all lattice points and the number of all cluster 
points respectively. 

Correlation lengths: £ is extracted as follows. 

Torus T 2 : Representing the point x as x = X\n\ + £2(^2 — n{) and using 
the vectors n x and n 2 from (|26| ) and (0), we actually compute 

C(x 2 )= £ C(x!,x 2 ). (53) 

xi=0 

In the range C 12 ^ ^, w e expect the following behaviour 

Xo ~~ 

C(x 2 ) = A cosh 2 2 , (54) 

where the constant A and the correlation length £ are extracted by fitting 
the Monte Carlo data. 

Sphere S 2 : In order to find the shape of the correlation function we start 
from a free, scalar, massive, Euclidean theory on the sphere S 2 . The propa- 
gator is as usual 

G — — - , (55) 
A + m 2 y ' 

A = ~d ll (Jg&'). (56) 

Using spherical coordinates, the following equation has to be solved: 

(57) 



or 



' -8 2 - -A-^(sin^) + (mR) 2 )G(<p, 0, <p r , & 



sin 2 1? ^ sin 1? 



(58) 

Lm 
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where the bar means complex conjugation. Thus, 

where $e[0,7r] denotes the angle between the points (<p,$) and ((p', i?'). We 
fitted the Monte Carlo data according to 



l-max Q7 _i_ 1 

C{0{x),O) = A^ -—zPticosti) (60) 

i=o 1(1 + 1) + (f ) 

in the range <C < 7T, truncating the series at / maa; = 100, which 

provides stable results. 



Scaling: Physical observables, e. g. Xt — (Q 2 ) /V an d £ > are measured in 
units of the average link length I, which is related to (3 through the 2 - loop 
renormalization group expression [H 



7? 



(61) 



_2 1 

Then, approaching the continuum limit Xtl and £1 are expected to behave 
like 



7?. 

47T „ \ 8tt^ 



(62) 
(63) 



(On T 2 replace I by the lattice spacing a .) 

We have checked the occurence of dislocations (i.e. lattice fields z a that 
have a topological charge Q = ±1 with respect to definition ([36]), (|37| ) but 

_2 

have an action below 87r/3/n i. e. the scaling exponent of xd ) by minimizing 
the action of the CP 1 model while keeping the topological charge fixed. We 
have found 





Sd 


sphere S 2 


2.09vr/? 


torus T 2 


2.06vr/? 



For higher values n the CP 1 configurations can be simply embedded. 
Therefore the dislocations found become suppressed for sufficiently large f3 
at n> 4 . 
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Let us consider the lattice CP 71 ^ 1 Model for large n on the torus T 2 . 
Then flBTD enables us to find the curve that £a _1 and Xt& 2 should lie on if we 
approach the infinite volume limit. Equations (|23|) and (|24j ) provide us with 



Xt (oo)a 2 = ±V« (aA L ) 2 (^*) . (64) 
n V Al J 

Introducing flBTD and 



= Sv/Sexp (-^(Invr + T'(l)) + ^ (65) 

into (|64|) gives x t (oo)a 2 . For n > 4 we expect our data to approach this 
curve for increasing n and £ 3> a. For the correlation length , we get 

= -W^'W-L-^. (66) 



40F aA L A Ai 



MS 



4. Results 

We have simulated the CP n ~ x model with n = 3,5, and 7, each case at a 
variety of (3 values and lattice sizes up to 10 2 and 60*120 , respectively, for 
T 2 and up to division = 6 for S 2 . Typically we carried out 10 4 cluster 
updates for thermalization followed by O(10 3 ) ( O(10 4 ) ) measurements of 
the topological charge (correlation function) separated by 10 ( 5 ) cluster 
updates. A cluster update means building up one cluster and flipping all the 
z's within this cluster. 

Figs. 5-9 show our main data. The errors of the topological susceptibility 
indicated were estimated by collecting the measured values into groups of 
0(100) single measurements and by taking the statistical error of the group 
averages. In order to determine the correlation length we splitted our data 
for the correlation function into 5 groups. For each of them we fitted the 
correlation length separatly and estimated the error from these values. 

In Figs. 5 we show the topological susceptibility for the cases n = 3, 5, 
and 7 on T 2 as a function of the " n -invariant" coupling / _1 = 47r(3/n . 
The approximate scaling behaviour or weak scaling violation seen for CP 2 
(Fig. 5a) is in accordance with the observations in and is related to the 
dominance of dislocations. In fact, at larger f^ 1 values the Monte Carlo 
time-histories for the topological charges have shown us the dominance of 
short-range fluctuations (i. e. only very short appearance of charges being 
unequal to zero). 
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For the cases CP 4 and CP 6 , where dislocations become definitely sup- 
pressed, we see a surprisingly strong deviation from asymptotic scaling in the 
opposite direction. This is in contrast to the conclusions drawn in [12] and 
in agreement with those of |14|]. The investigation in [12] concentrated on 
data produced in a smaller f~ l -range. Thus the falling tail at smaller bare 
coupling was missed. Nevertheless, the data presented there (with another 
lattice discretization) are compatible with ours. The time histories of topo- 
logical charges show freezing effects within sectors of non-vanishing Q over 
a few hundreds of consecutive charge measurements. This means long-range 
topological fluctuations occur. In spite of this, within the given (3 - or f^ 1 - 
range the cluster algorithm allows for sufficient tunneling between the sectors 
such that we consider our results for the topological susceptibility reliable 
up to the largest / _1 -values shown. 

In Figs. 5b, c, for comparison, we have drawn also the leading 1/n result 
for Xt(°°) according to (|4]). There seems to be a tendency to reach the 
1/n and continuum limit at very large (3 ! 

Figs. 6 present the analogous data for the case of the sphere S 2 . The 
observations and conclusions concerning the scaling violation are definitely 
the same. The 1/n limit cannot be given here because of the unknown 
A -ratio (compare with flBBD ). The comparison of Figs. 5 with Figs. 6 shows 
stronger finite-size effects in the S 2 -case than for the torus T 2 . 

In order to study the behaviour of the topological susceptibility as a 
function of the physical volume we measured the correlation length £ . Figs. 
7, 8 show the results of our fits to the correlation functions as described above 
and represented as functions of f~ l . Figs. 7a, b present the cases CP 4 
and CP 6 on T 2 , whereas Fig. 8 shows CP 6 on S 2 . All these data 
obtained for the correlation length are compatible with asymptotic scaling 
or even with slight scaling violations as reported in []nj. About possible 
convergence to the 1/n -limit we cannot draw here any conclusion. We 
need better statistics by at least one order of magnitude or/and an improved 
cluster, overrelaxation or multigrid algorithm Jl4[ . 



Finally, in Fig. 9 we study the volume dependence of the topological 
susceptibility for n = 7 and S 2 taking \t and the volume V (surface 
area of S 2 ) in units of the correlation length £ , which itself corresponds 
to the infinite- volume limit. The big error bars are due to the statistical 
errors of £ . Nevertheless, we see the expected tendency of a suppression of 
topological fluctuations in the limit of a small volume V . Fig. 9 contains 
the analytic result (|25"D , too, represented by a continuous line. Thus, we are 
qualitatively in the right ball park. But a real quantitative comparison with 
the continuum behaviour below V ~ £ 2 needs yet large efforts in future. 
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5. Conclusions 

In this paper we have studied the topological properties of the 2D CP n ~ x 
model on a lattice for varying n and for the cases of a (flat) torus and a 
sphere. We used the cluster update algorithm which can be economically 
implemented, and which allows to treat larger n 's and larger lattice sizes 
at intermediate coupling values than the standard Metropolis algorithm. We 
convinced ourselves that there is a reasonable tunneling rate between different 
topological sectors with long-living, large-scale field configurations in the 
coupling range considered. This does not mean that the algorithm allows to 
suppress critical slowing down. The contrary was shown in ||12|| . 

The main quantity considered was the topological susceptibility . It 
shows a strong violation effect from asymptotic scaling for n > 5 on T 2 
as well as on S 2 . In this case dislocations are known to become suppressed 
for sufficiently large (3 . On the other hand for n = 3 we recovered 
approximate scaling or slight scaling violation into the opposite direction 
due to the presence of dislocations. 

We compared the outcome of the Monte Carlo investigation with the 
analytical leading order 1/n result for infinite volume. Our data provide an 
indication that for very large f3 the 1/n result combined with asymptotic 
scaling could be reached. 

Unfortunately, within the given statistics the improved estimator method 
did not provide us correlation lengths £ for large volumina with sufficient 
high accuracy in order to study the limit of small volume V — > in terms 
of physical units £ . So the comparison of the semi-classical and 1/n 
behaviour on the sphere S 2 in this limit deserves further investigations. 
So far we obtained qualitatively a reasonable behaviour at V > £ 2 . The 
multigrid method can provide a possibility to improve that. We hope to 
come back to this question very soon. 
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Figure captions 



Fig. 1: First three steps of refinement for discretizing the sphere S 2 . 

Fig. 2: Link length distribution on S 2 for refinement step N = 5 . 
Average link length is I = .08794 for a unit sphere. 

Fig. 3: Distribution of cluster sizes produced with the one-cluster algorithm 
[2^ for T 2 with lattice size 100 x 100 and CP 6 at (5 = 4.6 . The average 



size is 6516. 

Fig. 4: Comparison of numerical < Sij > values with analytic ones 
(continuous lines) at weak coupling. Squares (triangles) correspond to CP 2 
(CP 4 ). 

Fig. 5: Topological susceptibility as a function of f~ l for the torus T 2 
and CP 2 (a), CP 4 (b) and CP 6 (c), respectively. The straight line 
corresponds to an asymptotic scaling behaviour. The dashed line is the 1/n 



result according to (]64|) 



Fig. 6: Topological susceptibility as a function of / _1 for the sphere S 2 
and CP 2 (a), CP A (b) and CP 6 (c), respectively. The straight line 
corresponds to asymptotic scaling. 

Fig. 7: Correlation length £ as a function of / _1 for T 2 and CP 4 (a) 
and CP 6 (b), respectively. Straight and dashed lines as in Figs. 5. 

Fig. 8: Correlation length £ as a function of / _1 for S 2 and CP 6 . 

Fig. 9: Topological susceptibility as a function of the physical volume, both 
expressed in units of the correlation length. The continuous line shows the 



one-instanton result ( p5|) acc. to [16|. The dashed line indicates the infinite 
volume 1/n result, cf. Eq. (p3|). 
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